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Abstract. A DPMJET-III model (DPMJET-III-2011) with chain fusion adjusted 
to include energy. dependent parameters is used to calculate inclusive distributions in 
p-p collisions at LHC energies. Presented are charged hadrons rapidity distributions, 
transverse momentum distributions, multiplicity distributions as well as multiplicities 
at mid-rapidity as function of the collision energy. For hadrons with strangeness we 
present cms-rapidity distributions and transverse momentum distributions. With the 
considered merely energy-dependent adjustments the obtained agreement with the 
transversal A and S distribution is not satisfactory. 
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1. Introduction 

Monte Carlo codes based on the two-component Dual Parton Model involving soft 
and hard hadronic collisions producing chains of particles are available since almost 20 
years pTj. The present codes are: 

PHOJET for h-h and 7-h collisions [2] 
DPMJET-III for h-h, h-A and A-A collisions |3] 

In distinction to earlier versions DPMJET-III is based on PHOJET for its h-h 
collisions. In such collision it is therefore - except for a few additions like the fusion 
discussed below - identical to PHOJET. PHOJET describes the production of strings. 
For the string decay it calls PYTHIA version 6.412 For a few special cases we found 
it necessary to change the PYTHIA fragmentation. These were done in the DPMJET 
part, leaving the PYTHIA code itself untouched. 

As we use the full program we will refer below just to DPMJET-III. We now outline 
its main additions. 

Comparing DPMJET-III to RHIC data it was learned that something had to 
be done to decrease the particle density. As the strings are quite dense in impact 
parameter space interactions between strings are plausible. The expected percolation 
was modeled as fusion of close hadronic chains implemented in DPMJET-III [5] in 2004. 
The obtained reduction was very essential for central collisions of heavy ions, but fusion 
also changes the particle production in very high energy p-p collisions when the number 
of contributing chains obtained by a Glauber / eikonal formalism gets sizable. 

RHIC and Fermilab data also contain interesting information about particle 
antiparticle ratio's [6j. For the baryon/antibaryon distribution the string fusion 
mentioned above can be significant (p.e. two quark-antiquark strings can fuse to a 
diquark-antidiquark string yielding baryons and antibaryons). 

In the diquark string decay used in PYTHIA one observes a dip in the ratio of 
the fl/Cl spectra not seen in the data. A solution of the problem is to include a small 
contribution of diquark-diantiquark mesons production in the first rank so that Q can 
appear in the second rank. The idea is that such tetra-quark mesons are always produced 
but decay too fast to be identified in mass plots. 

The LHC experiments did compare DPMJET-III to particle production at LHC- 
energies, see [7], [8] and [9]. There were some successful predictions of DPMJET. 
However, LHC experiments found that around 7 TeV the multiplicity rises faster with 
energy than predicted by DPMJET-III. 

In order to make the program usable for ongoing data analyses at LHC energies 
we adjusted the program to improve the agreement with available experimental results. 
We allowed for an energy-dependence of string decay parameters. No differentiation 
between softer and harder strings was attempted. The new results of this modified 
version will be reported in section 3 and 4. 
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2. Modifications of DPMJET-III needed for LHC energies 

There are essentially three additional modifications of DPMJET-III implemented in 
order to get better agreement with LHC data on particle production in p-p collisions. 

(1) The first modification is connected to a problem with collision scaling known since 
2004 [TU] . DPMJET-III uses an eikonal formalism to determine the size of various 
multiple scattering contributions ^^ty where n is the number of chains and 

a(^^ is the Regge intercept depending on the diquark, valence quark or sea quark 
nature of the forward/backward parton of i-th chain. Let us consider the forward 
direction. For each such configuration the attributed energy fractions {xi} to these 
partons are then chosen with a factorizing structure function of the form: 



in which the energy available for a scattering process depends on the reminder. 
There are two kinds of chains in DPMJET. Hard chains produced by hard 
collisions of partons from the colliding hadrons (typically large p±) and soft chains 
representing soft hadron production in the collisions. In the factorizing formalism 
soft processes affect the energy sampled in hard processes. 

This turned out to be too simple. Experiments [H] gave evidence for collision 
scaling in not too central scattering processeiji- Collision scaling means, that exactly 
as many hard chains are produced as predicted by considering just hard collisions. 
To correct for the missing collision scaling an additional parameter was 
introduced [lO] which increases the number of hard collisions in such a way that 
collision scaling is obtained. We here adjusted these constants to the LHC data. 

(2) LHC experiments [3 [3 [9] found that the produced hadron multiplicity rises 
with energy somewhat faster than the model. To obtain this faster rise with 
collision energy one needs to introduce energy- dependent parameters. This is not 
unreasonable. The time scale of the initial scattering is inversely proportional to 
the energy. This causes a more localized string and a widening of the distribution 
of the string ends which was observed as on contribution to the multiplicity 
dependence of the average transverse momentum < p± >„. There are of course 
many different parameter in PYTHIA which can be tuned in an energy-dependent 
way. 

The solution which was adopted for the moment is not to search for a suitably, 
possibly theoretically plausible energy-dependence but just adjust parameters. The 
two parameters which determine the multiplicity of fragmenting chains are the Lund 
parameters PARJ(4i) and PARJ(42) for which we use the following values: 




I For central heavy ion collision collision scaling is lost as transport effects become important. So far 
these effects are not implemented. 
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In ([T]) and ([2]) we do not continue to change PARJ(41) and PARJ(42) for i?cm larger 
than the maximum LHC energy of 14000 GeV. For the moment (in "DPMJET-III- 
2011") we replace in the FORTRAN code PARJ(41) and PARJ(42) by the values 
given in ([1]) and ([2]) as soon as in the input cards a change of PARJ(41) and 
PARJ(42) is demanded. 

(3) The third modification is connected to the production of strange hadrons. The 
production of mesons and of A and S~ hyperons in p-p collisions was measured 
by the CMS Collaboration [19]. The program gave more than measured by CMS, 
while it obtained less A and S~ production than measured by CMS. To increase the 
agreement with the measurements in this regard more energy-dependent parameters 
have to be introduced. 

Hyperon and strange meson production in DPMJET-III is controlled by the 
Lund parameters PARJ(l), PARJ(2), PARJ(3), PARJ(5) and PARJ(6). The 
default of the parameter PARJ(2) was not touched. For the other parameters 
the following energy-dependent values (in the energy range up to Ecm = 14 TeV) 
were implemented: 
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3. Comparison of DPMJET-III-2011 results with LHC Data on charged 
hadron production 

We start to discuss the nsd (non single diffractive) and inel (inelastic) pseudo-rapidity 
distribution dNch/drjcm measured by the CMS and ALICE Collaborations. 



dndetcms2orig2parj1 3 



11 

10 - 
9 - 
8 - 

7 - 

.J 

6 - 

J 

5 - 
4 - 



11-2011 nsd 14000 


GeV ■ 


inel 


14000 


GeV ■ 


nsd 


7000 


GeV ■ _ 


inel 


7000 


GeV ■ 


nsd 


900 


GeV ---M--- 


inel 


900 


GeV ■ 


CMS nsd 


7000 


GeV ; X- -; - 


nsd 


900 


GeV ; ^ ; 


ALICE nsd 


900 


GeV + 


inel 


900 


GeV ^ 
























. - - • 


■■•■.„. 







IX- 



X. 



2 I I I I I I I 

-3-2-10123 

Fig.l Central rjcm distributions of charged particles in \fls) = 900, 7000 and 14000 
GeV p-p collisions compared to non single diffractive (nsd) and inelastic (inel) pseudo- 
rapidity distributions obtained with DPMJET-III-2011. The experimental data for nsd 
and inel collisions are from the CMS (nsd) [12] and from the ALICE [13] collaboration. 
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In Fig.l and Fig. 2 we present for p-p collisions the non single diffractive data 
from CMS [12] at 900, 2360 and 7000 GeV and the non single diffractive and inelastic 
data from ALICE p3| at 900 and 2360 GeV and compare them with the results from 
DPMJET-III-2011. Excellent agreement is obtained. 

Also included are the results at 9000 and 14000 GeV in p-p collisions. At these 
energies the distributions are expected to be be measured at the LHC in future. 



E 6 

o 



o 

-z. 

-o 5 



dndetcms2orig2parj1 1 



8 - 



7 - 



4 - 



3 - 



DPMJET-III-2011 nsd 2360 GeV 

9000 GeV 
inel 2360 GeV 
9000 GeV 
CMS nsd 2360 GeV 
ALICE nsd 2360 GeV 
inel 2360 GeV 



+ 



+ 



+ 



+ 



+ + 



" + 







Fig. 2 Additional central r/cm distributions of charged particles in y (s) = 2360 and 
9000 GeV p-p colhsions compared to non single diffractive (nsd) and inelastic {inet) 
distributions obtained with DPMJET-III-2011. The experimental data for nsd and inel 
collisions are from the CMS insd) [12j and from the ALICE [13] collaboration. 
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The energy-dependence of the central density dN/drjcm at rjcm = is presented 
in Fig.3 for p-p coUisions of nsd and inel events. The DPMJET-III-2011 resuhs are 
compared with data from various energies. In all cases a good agreement is obtained. 
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Fig.3 Central rjcm values of all charged particles in p-p collisions compared to 
DPMJET-III-2011 results. The experimental data are from the CMS collaboration [12] 
for nsd collisions and from the ALICE collaboration [13] for nsd and inel collisions. 
Further data are from UA5 [14], the ISR [15], STAR [I6], UAl [IT] and CDF [H]. 

In Fig.4 we compare pt distributions from the DPMJET-III-2011 in p-p collisions 
at ^/s = 900, 2360 and 7000 GeV with representative experimental data points from the 
CMS Collaboration [7]. The agreement between the modified program and the CMS 
data points is good. 
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Fig. 4 Transverse momentum distributions in p-p collisions at ^/s = 900, 2360 and 
7000 GeV. We compare experimental data from the CMS Collaboration [7] to DPMJET- 
III-2011 calculations. 

In Figs. 5-7 we compare the multiplicity distributions for |?7| < 1. with experimental 
data from the ALICE Collaboration [13] . Again a reasonable agreement is obtained. 
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Fig. 5 Multiplicity distributions in p-p collisions at ^/s = 7000 GeV. We compare 
data for Ir^l < 1. from the ALICE Collaboration [13] to the DPMJET-III-2011 result. 
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Fig. 6 Multiplicity distributions in p-p collisions at ^/s = 2360 GeV. We compare 
data for \r]\ < 1. from the ALICE Collaboration [13j to the DPMJET-in-2011 result. 
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Fig. 7 Multiplicity distributions in p-p collisions at ^/s = 900 GeV. We compare 
data for Ir^l < 1. from the ALICE Collaboration [13] to the DPMJET-III-2011 result. 



4. Comparison of DPMJET-III-2011 results with LHC data on the 
production of strange hadrons 

Strange hadron production in p-p collisions was measured by the CMS Collaboration 
[T9] . It determined the production of mesons and the production of A and S~ 
hyperons. 

Similar data on the production of strange hadrons were also given by the ALICE 
Corporation [20j . We did not include them so far as they do not affect the consideration 
discussed below. 

In Fig. 8 we compare the production of and A hyperons in dependence of the 
dn/dycm ■ With energy-dependent parameters good agreement of the DPMJET-III-2011 
and the CMS measurements is obtained. 

But the situation is not perfect. CMS also measures transverse momentum 
distributions. Comparing transverse momentum distributions we find the shape of the 
distributions to differ. This can be seen in Fig.9 showing the transverse momentum 
distributions of A hyperons. Above 1 GeV the model is below the data. A similar 
problem seems to appear in many model calculations pT] . 

If we would also adjust the parameters in such a way, that the agreement between 
the transverse momentum distributions is optimal in this region, we would obtain a 
disagreement in the dn/dycm distributions. 
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Fig. 8 dn/ducm distributions in p-p collisions of and as well as A and A. We 
compare LHC data from the CMS Collaboration |19j to DPMJET-III-2011 result. 
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Fig. 9 Transverse momentum distributions of A and A. We compare experimental 
data from the CMS Collaboration [H] to DPMJET-III-2011 result. 

Unfortunately here the situation becomes even more problematic. In Fig. 10 we 
compare the dn/di/cm distributions of S hyperons in the DPMJET-III-2011 with the 
measurements of CMS. The modified model predicts S distributions about three times 
as large as measured by CMS. We have modified the parameters in such a way, that 
the A hyperons agree with the CMS data. The same parameters should also lead to 
agreement for the S hyperons. They do not. We can only conclude, that so far we do 
not fully understand the production of H hyperons in DPMJET-III. 
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Fig. 10 dn/dycm distributions of S and S. We compare experimental data from the 
CMS Collaboration [19j to DPMJET-III-2011 calculations. 



5. Conclusions 



DPM JET-111 is a code for hadron production in hadron-hadron, photon-hadron, hadron- 
nucleus and nucleus-nucleus collisions [21 |3] which is about 10 to 15 years old. The 
measurements of hadron production in p-p collisions at LHC energies gave the occasion 
to check how well DPMJET-III agrees to the data in this higher energy region. We 
knew already from comparisons at the lower FERMILAB energies, that not all features 
of DPMJET-III are valid at higher energies, a well known example is the collision scaling, 
which was found at FERMILAB in hadron-nucleus and nucleus-nucleus collisions and 
which is not a property of DPMJET-III. 

Comparing DPMJET-III with the LHC data we found further problems of 
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DPMJET-III. The energy- dependence of hadron production measured by the LHC 
collaborations at 7 TeV did differ from the one predicted by the original DPMJET-III. 
In the present paper we find, that this energy-dependence can be corrected in DPMJET- 
III by making in DPMJET-III the PYTHIA-parameters energy dependent. We consider 
the introduction of energy-dependent parameters only as a temporary solution to get 
agreement of DPMJET-III with the new LHC data. 

A more permanent solution will require deeper changes in the program. Of course 
the parameters entering DPMJET-III, PHOJET and PYTHIA should not be adjusted 
per hand to determine the energy- dependence of the hadron production models. This 
energy- dependence should be an intrinsic property of the hadron production models. 
We conclude that we need a new version of the model which agrees better with the data 
in the new energy region opened by the LHC. 
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